Simple pulses for elimination of leakage in weakly nonlinear qubits 
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In realizations of quantum computing, a two-level system (qubit) is often singled out from the 
many levels of an anharmonic oscillator. In these cases, simple qubit control fails on short time 
scales because of coupling to leakage levels. We provide an easy to implement analytic formula 
that inhibits this leakage from any single-control analog or pixelated pulse. It is based on adding 
a second control that is proportional to the time-derivative of the first. For realistic parameters 
of superconducting qubits, this strategy reduces the error by an order of magnitude relative to the 
state of the art, all based on smooth and feasible pulse shapes. These results show that even weak 
anharmonicity is sufficient and in general not a limiting factor for implementing quantum gates. 
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Quantum information processing devices are paving 
new inroads into the understanding of coherent processes 
and their applicability to high-complexity computation. 
The fundamental building blocks for these devices are 
quantum bits (qubits), which are quantum two-level sys- 
tems. Naturally occurring qubits such as spin g particles 
and photons are difficult to both isolate and control. For- 
tunately, effective qubits can also be manufactured from 
an anharmonic multi-level system. These systems can be 
understood as a generalized atom, and all share the com- 
mon problem of leakage out of the qubit subspace if the 
control bandwidth is comparable to the anharmonicity. 
Examples of such systems are superconducting qubits [1- 
15], optical lattices [16], and quantum dots [17]. 

In superconducting qubits [1] , researchers have demon- 
strated single qubit gates [10, 11], two qubit gates [4, 5], 
and the use of a quantum bus [6, 7], suggesting the possi- 
bility of a scalable technology. Decoherence has dropped 
dramatically [9, 14, 15] from the original designs [2], due 
to both device and operation improvement [3, 9, 12]. To 
date, the single-qubit accuracy of such systems has been 
in the range of 98% [10]- 99.3% [11], which is a major 
achievement, yet not sufficient for quantum computing. 

For large scale computation, gate errors are required 
below a fault-tolerance threshold, conservatively esti- 
mated at 10 -4 [18]. In superconducting qubits such as 
the Phase [8, 9, 10] and the Transmon [11, 12, 13, 14] 
qubits, long pulses are limited by decoherence and short 
pulses by excitations out of the qubit space [11, 19]. 
Thus, one prevalent approach to achieve high fidelity 
gates is to increase the coherence time; alternatively, the 
approach we will take is to reduce the gate time by using 
optimized control pulses that avoid leakage. 

Thus far, trying to minimize gate times through bet- 
ter control has been met with limited success. Little im- 
provement has been achieved beyond what is possible 
with a Gaussian waveform [19, 20]. Optimal control the- 



ory with only amplitude modulation has required pulses 
which have sharp features [21, 22, 23]. For example, the 
pulses in [19, 21] require 3 pulse steps and 2 waiting pe- 
riods, which requires currently impractical pulse shaping 
hardware. Furthermore, these pulses are limited to a 
minimal time, given by the inverse of the anharmonicity. 

We propose to solve this problem using two quadra- 
ture controls for effective phase modulation, which al- 
lows us to achieve instantaneous cancellation of leakage 
using slowly varying controls. We present an analytical 
approach to accomplish this called Derivative Removal 
by Adiabatic Gate (DRAG). Then, we map the full po- 
tential of the technique using numerical refinement and 
consider the case when decoherence is included. 

The system. — We consider the lowest three levels of 
a driven slightly anharmonic energy spectrum with only 
nearest level coupling; the first two levels comprise the 
qubit and the third level accounts for leakage, provided 
its population is minimized during the entire pulse. That 
is, in the lab frame, the Hamiltonian can be written as 



H = h 



(1) 



where IF,- = \j)(j\ is the projector for the j level, 
aj = \j — is the lowering operator, and the cor- 
responding transition energies are denoted hu>j , with the 
ground energy set to zero. The intrinsic anharmonicity 
of the system is A = ll>2 — 2cji . The parameter Xj weighs 
the relative strength of the — 1 to 1 — 2 transition, and 
without loss of generality we will take Ai = 1 and X2 = X. 
Driving and control of the system is represented by 



£{t) = 



£ x (t) cos(Lu d t) + £y (t) sin(cj d t) , < t < t g 
0, otherwise, 



(2) 

which is a single frequency carrier with two independent 
quadrature controls, £ x (t) and £ y (t); t g is the time taken 



2 



for one gate operation. 

Making the rotating wave approximation (RWA) and 
moving to the interaction frame with respect to the drive 
frequency uj^, the Hamiltonian can be written as 



H 1 



j=l,2 



£ x (t) 



V7- 



£v(t) 



where a x k 



Si 



|fc)(iH 

Wi — and ^2 



(3) 

and<7j !fe =<|*)01-i|j)(fc|- 
= A + 26i are the detunings of the 



transitions with respect to the drive frequency u>d- When 
the drive frequency is resonant with the qubit frequency 
(Si = 0), S2 = A is the anharmonicity. 

In the ideal case, A = 0, the Hamiltonian and its evo- 
lution reduce to that of a qubit. Thus qubit operations, 



t/ ideal = Texp -- / H R dt =e*(C/ qb + e^n 2 ), 

(4) 

can be trivially driven by the controls, Eq. (2), since 
the effect is t/qb, which is a unitary that acts only in the 
qubit subspace, and the global phases <pj are irrelevant 
[21]. For example, the NOT gate (which without loss 
of generality we choose to focus on) is implemented by 
simply setting £ x = £ Vi £ v = 0. Then, the time-ordering 
operator T can be dropped and any controls that satisfy 
J a £ n dt = 7r lead to the desired gate t/ q b = Cq x . 

However, in the more general case when A is non-zero, 
the pulse can result in large unwanted leakage into the 
third level. Fourier analysis can show how this leakage is 
dictated by the anharmonicity A [19, 24]. The standard 
approach to diminish this error is to use Gaussian (£g) 
[20] or tangential (£t) envelope shaping 



£ G {t) = j^exp 

£ T (t)=U 



{t-t g /2f 
2a 2 



tanh I — I — tanh ' 

a I \ a 



D 

to -I 



D 



(5) 



where a is the standard deviation for £q and the rise- 
time parameter for £t- A is chosen such that the correct 
amount of rotation is implemented (e.g. tt for a NOT) 
while B enforces the pulse to start and end at zero. The 
small frequency response bandwidth of the Gaussian en- 
velope ensures little excitation at the leakage transition 
frequency. Additionally, Gaussians and tangential pulses 
have relatively simple pulse-shaping hardware require- 
ments. This being said, these pulses still suffer large 
errors at short gate times, as shown by the solid blue 
line in Fig. 1, where the error is 1% at a gate time of 
6 ns (vertical dashed line), and scales as (£,n ax ) 2 / A. To 
quantify this error, we use the gate fidelity averaged over 
all input states existing in the qubit Hilbert space, which 
using an argument similar to [25] gives 



F„ 



j=±x,±y,±z 



Tr(Ui de ,i Pj U^Mp( P] ))) (6) 



where Mp(pj) is the actual process in the three dimcn- 
sionional Hilbert space, pj are the six axial states on the 
Bloch sphere, and £/idcai is defined in Eq. (4). 
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FIG. 1: Error vs. pulse duration for Gaussian (a=0.5t g ) (blue 
solid), Gaussian (a—0.5t g ) with DRAG (red dot-dashed), 
tangential (a=0.6t s -1.4) with DRAG (green dashed), and 
GRAPE pulse (purple crosses with dotted guiding line). The 
standard Gaussian and DRAG pulses are continuous, while 
the GRAPE pulse has a pixel size of Ins. Parameters are 
A = v 7 ^, A = 2tt(-400 MHz) (typical for the Phase [8, 9, 10] 
and Transmon [11, 12, 13, 14] qubits). The top axis is set to 
dimensionless units, 1/|A|. 

Derivative Removal by Adiabatic Gate (DRAG). — To 
go beyond these simple pulses we introduce an adiabatic 
transformation V that allows us to work entirely in the 
qubit subspace. This transformation is 



V(t) = cxp [-i£* + Xal 2 ) /2A] , 



(7) 



where £ x / A is the adiabaticity parameter. For the gates 
to be identical in both frames we require 1^(0) = V(t g ) = 
1, such that by the end of the pulse it has no net effect. 
This is achieved by setting £ x (0) = £ x (t g ) = 0. In this 
frame, to first order in £ x /A, the effective Hamiltonian 
H v = VH R Vi + ihVV^ reads 



H v /h: 



£ J 



(A 2 -4)£2 
(*i )Hi 



\C2 

x _x 1 fx 

r £y i- 



(A 2 + 2)£\ 



)n 2 



2A 



4A 

(°o,i + ^1,2) 



(8) 



Here, we see we can remove the imaginary incrtial term 
and the ac-Stark shift (phase) error with the second 
quadrature and drive detuning set to 



mi 



(9) 



Thus, we eliminate the leakage (second line of Eq. (8)) 
to order £^/A 3 , which in this frame comes about from 
the induced — 2 transition [28]. This error can also 
be removed by adding the term — i\£ 2 &Q 2 /8A 2 to the 
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transformation V. This results in a higher order error 
being added to the — 1 transition, 3££<7q 1 /16A 2 , which 
can easily be compensated by rescaling £ x . Increasing the 
precision of the transformation up to 5th order results in 



£ x {t) 



(A 2 -4)£3 (13A 4 -76A 2 + 112)£| 



8A 2 



128A 4 



cy/,\ . 33(A 2 — 2)5 
£ {t)= ~ A + 24A^ ■ 



5 1 (t) 



(A 2 - 4)£ 2 (A 4 - 7A 2 + 12)£ 4 



4A 



16A 3 



(10) 



These pulses remove most of thc leakage error, as 
demonstrated in Fig. 1, often with only the first order 
of corrections. The green dashed (optimized tangential) 
and red dot-dashed (optimized Gaussian) lines in Fig. 
1 show that DRAG reduces the gate error by orders of 
magnitude (e.g. 4 orders at 6ns). In fact, we verified with 
randomized testing that DRAG is largely impervious to 
envelope shape, provided it has the right area and starts 
and ends at 0. Fig. 2 A shows an optimized Gaussian. 

If Si cannot be controlled in real time, phase ramping 
can be applied [26]. It amounts to transforming to the 
frame given by transformation P = exp[— i6(t) ^2 j^-j] 
where 6(t) = Jq(Si(s) — 5~i)ds . Then, alternatively, we 
can pick 6{ = 8± = f f* a 5i{t)dt, £ x (t)' = £ x {t) cos 0(f)- 
£^(f) sin 0(f), and £ v {t)' = £ y (t) cos 0(f) + £ x (t) sin 0(f). 
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FIG. 2: A. 6ns analog Gaussian (a = t g /2) with DRAG; 
B. 4ns digital GRAPE with no detuning and pixels at Ins. 
In both cases, the blue solid line is the £ x control and the 
green dashed is the £ y control. The red dot-dashed line in 
A represents explicit detuning, Si. Other parameters are the 
same as in Fig. 1. 

Numerical optimization ( GRAPE). — We apply opti- 
mal control theory to our system by employing the gra- 
dient ascent numerical optimization algorithm, GRAPE 
[27]. Wc try to minimize the gate error, 1 — F g , with F g as 
defined in Eq. (6) . This is equivalent (for dissipation-free 
evolution) to maximizing 



TtKOl^dcaltWlO) + (ll^C/g.apell}! 2 ), (11) 



where U glape is the unitary evolution found by GRAPE 
[21]. GRAPE works by time-slicing the evolution into 



discrete time pixels of constant amplitudes £ x and £ v , 
as in Fig. 2B. The gradient of $2 with respect to each 
control at each time step provides a direction in search 
space. That is, each iteration of thc algorithm increments 
the controls at each time step of the pulse proportionally 
to the gradient with respect to the control at that time 
(including the sign). This unique search direction allows 
for fast optimization, particularly with few pixels, and 
the algorithm quickly converges to an optimal solution 
(1 - $2 < 10~ 5 ) for all gate times (see Fig. 3). 

Restricting GRAPE to Ins pixels (which is consistent 
with current experimental limitations [29]) and setting 
the initial condition for thc £ x quadrature to be the 
Gaussian pulse, Eq. (5), the algorithm quickly converges 
to the optimal solution. After optimization, the shape of 
the Gaussian in the £ x quadrature is largely unchanged. 
However the £ y quadrature changes from to a shape 
similar to thc derivative of the Gaussian, sec Fig. 2B. 
GRAPE essentially takes the idea of DRAG to infinite 
order and discretizes the pulse. The shapes of the pulses 
are smooth with no sudden rises or falls in the control 
amplitudes. Thc purple dotted line in Fig. 1 plots the 
gate error as a function of gate time for GRAPE. The ex- 
tra optimization increases the performance over DRAG. 

The two quadrature control GRAPE pulse outper- 
forms its one quadrature control version [21]. Fig. 3 
shows the minimal pulse duration to obtain gate errors 
smaller than 10~ 5 vs the pixel width for both one (red 
solid) and two (blue dashed) controls. Here, we see that 
at small pixel widths (continuous limit) the one control 
saturates at a time of 2tt/A [21] (black dashed line), 
whereas the two controls allows for arbitrarily small pulse 
times. This is because, in the one control case, phase can- 
celation of the leakage is done through phase accumula- 
tion by the natural precession of thc third level, whereas, 
for two controls, it can be done instantaneously. The in- 
sert of Fig. 3 shows such a continuous two control pulse; 
this pulse clearly displays a combination of DRAG and 
composite features [30] . At the larger pixels (more exper- 
imentally realizable), the two controls offer a substantial 
improvement over one control. For example, at Ins pix- 
els, there is a factor of 3.5 improvement. 

Decoherence. — In realistic physical systems, decohcr- 
ence cannot be neglected. Here, we include this effect by 
simulating the evolution by a master equation for each 
optimal pulse. The master equation we use is 



.7=1,2 



XV[aj)p + X ViU^p 



(12) 



where p is the density matrix, T( and represent re- 
laxation and pure dephasing times respectively, and V is 
the damping superoperator, defined as 

V\A\ P = ApA^ - \tfAp - \pA^A. (13) 
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FIG. 3: Pulse duration vs. pixel width for 1 control (blue 
solid) and 2 controls (red dashed). Pulse duration is calcu- 
lated as the minimal time to bring the error down to 10 -5 . 
The insert shows the two control GRAPE result for t g — 1/4 A 
(blue solid is £ x , and green dashed is £ y ). Other parameters 
are the same as in Fig. 1. 

For the simulation, wc used parameters consistent with 
the Transmon qubit and as such take pure dephasing to 
be zero. The results are shown in Fig. 4. Fig. 4A plots 
the gate error (1 - f s ) as a function of gate time for 
T\ = 40/iS. Here, we see there is an optimal value for the 
gate time (when the error due to decoherence is about the 
same as the error due to leakage) and for both the DRAG 
(green dashed) and the GRAPE (purple dotted) pulse, 
this is much less then the standard Gaussian pulse (blue 
solid) with an optimal error of 5.3 x 1CK 5 and 3.4 x 10~ 5 
respectively. Fig. 4B plots this minimum as a function of 
T\ where we observe that at current experimental values 
Ti « (1 -4) fis [14] DRAG and GRAPE still out preform 
the simple Gaussian by about an order of magnitude. 

Conclusions. — The analytical and numerical results 
presented here demonstrate that it is possible to elimi- 
nate leakage in physical qubits by using a second quadra- 
ture control approximately equal to the derivative of the 
first. The predicted single-qubit error rates are extremely 
low, even when we include decoherence. While in this pa- 
per this procedure is only applied to leakage, it will also 
find application to other systems, for example turning 
off unwanted coupling in multi-qubit and qubit-oscillator 
systems for single and two qubit operations. 
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